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We investigate theoretically the efficiency of intra-molecular hot carrier induced impact ionization 
and excitation processes in carbon nanotubes. The electron confinement and reduced screening lead 
to drastically enhanced excitation efficiencies over those in bulk materials. Strong excitonic coupling 
favors neutral excitations over ionization, while the impact mechanism populates a different set of 
states than that produced by photoexcitation. The excitation rate is strongly affected by optical 
phonon excitation and a simple scaling of the rate with the field strength and optical phonon 
temperature is obtained. 



The excellent electrical properties and direct gap of 
carbon nanotubes (CNTs) offer the possibility of a uni- 
fied electronic and optoelectronic technology based on 
this material . The electroluminescence (EL) in CNTs 
field effect transistors has already been demonstrated, 
where electrons and holes were independently injected in 
an ambipolar device |^, Recently, light emission 
has been observed under unipolar transport conditions 
in the suspended tubes and from the CNT-metal con- 
tacts 0, 01 . In the suspended sample the light intensity 
is stronger by a factor of ~1000 compared to ambipo- 
lar devices. The mechanism of this enhanced EL was 
proposed to involve an unusually efficient intra-nanotube 
impact excitation process involving the hot carriers in 
CNTs, generated by the high local fields at the interface 
between the suspended and nonsuspended parts of the 
tube and the metal-nanotube contact. By utilizing such 
impact excitation process in quasi one-dimensional con- 
fined structures, high exciton densities can be achieved 
for probing the interactions of such boson systems and 
for generating ultra bright nanoscale light sources used 
in nanophotonics. 

Here, we provide detailed theoretical study of the im- 
pact excitation in semiconducting CNTs. In low dimen- 
sional materials, the Coulomb interactions are weakly 
screened, and large exciton binding energies in CNTs 
have been predicted theoretically 0) S U and verified 
experimentally Ull . This presents an additional 
challenge for the computation of impact excitation rates 
in CNTs. We find that the impact scattering rate is very 
fast, up to five orders of magnitude larger than those in 
bulk semiconductors ^3 ■ This can be easily rationalized 
since the same Coulomb interactions responsible for the 
large exciton binding are also responsible for the impact 
excitation scattering. The electron-hole correlations in 
the final state make the impact excitation process, on 
the average, twice more efficient than the impact ioniza- 
tion, where the created electron hole pairs are treated 
as independent particles. Impact excitation generates a 
different initial distribution of excited states than that 
produced by photon excitation due to the different selec- 
tion rules and momentum conservation requirements. 

It has recently been shown that under high bias the 



energetic carriers in nanotubes excite optical and zone- 
boundary phonons 0, 0, 0, 0| and generate a non- 
equilibrium phonon distribution 0, 0], particularly 
when energy dissipation to the substrate is suppressed 
as in suspended CNTs We have calculated the ef- 

fect of optical phonon excitation and found that the ex- 
citon production rate depends exponentially both on the 
hot phonon temperature and on the electric field, and 
the effect can be described by a simple scaling relation. 
As a result, optical phonon excitation significantly in- 
creases the EL intensity under unipolar transport con- 
ditions. The energy spectra of the excitons generated 
consist of several peaks corresponding to excitons with 
different angular momenta L: Efi, E12, and E22, where 
E^j stands for exciton with electron and hole primarily 
from bands i and j. The relative intensities and the over- 
all excitation efficiencies are functions of the electric field, 
the optical phonon temperature, and the CNT chirality. 

The impact ionization scattering rate of an electron 
with momentum and energy (fci, e^J scattered to state 
{ki — q, £%-^_q) plus an e-h pair with a hole (— fc2, -e^^) 
and an electron {k2 + q, , ) is given by |20j : 
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where i' = kik2q, — J{kik2, k2 + qki — q) and — 
J{kik2,ki — qk2 + q) are the direct and the exchange 
integrals: 

J(fcifc2, fci -qk2 + q) = J dfidr2KkA''i)'^*vk2ir2) 
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V{fi,f2) is the Coulomb potential screened by the di- 
electric environment of the nanotube and ^'^(r) and 
^'*(r) are the electron and the hole wavefunctions, given 
by a TT-orbital tight-binding model with hopping matrix 
element t = 3.0 eV. In Eq. ^ we neglect e-h pair oc- 
cupation in the initial state and, hence, processes where 
electron can absorb an e-h pair. 

To include the effect of e-h correlation in the final 
state we follow I2H and show that the impact excita- 
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tion scattering rates due to the production of singlet (S) 
and triplet (T) excitons are given by: 
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^^^9. = E^^^.^^^'c., (3) 

fc2 

where E^^ and are the singlet and triplet exciton 
transition energies 22|, is the solution of the Bethe- 
Salpeter equation 8, 9] for the exciton two-particle wave- 
function: 

^^qiri,r2) ^J2'^kq^ck+q{ri)nkir2). (4) 
fc 

We model the screening of the Coulomb interaction by 
the dielectric constant of the medium embedding the nan- 
otube ,9|, both for the Bethe-Salpeter equation kernel 
and the impact ionization Coulomb potential Eq. ||2J). 
The CNT diameter dependence of exciton binding ener- 
gies obtained by two-photon florescence excitation spec- 
troscopy [23 agree very well with our model calculations 
by choosing e = 3.3, a value which we use in the rest 
of the paper. 

The impact ionization and impact excitation scatter- 
ing rates for a tube of diameter d = 2 nm are shown in 
Fig.^ and^ respectively. This is a typical diameter of 
tubes used in opto-electronic applications and grown by 
chemical vapor deposition. The impact excitation rates 
in CNTs are roughly twice as large as the impact ioniza- 
tion rates due to the e-h confinement in the final state. To 
understand the structure of the impact ionization scatter- 
ing rates in Fig. we need to consider the implications 
of energy and momentum conservation. The momen- 
tum has two components: a discrete angular momentum, 
which labels the bands, and a linear momentum along the 
tube axis. The CNT bands are non-parabolic and an an- 
gular momentum conservation law plays a crucial role in 
determining the threshold energy of the impact excitation 
|24l |. As seen from Fig.^, angular momentum conserva- 
tion allows electrons in the third, 63, and the forth, 64, 
bands to undergo impact ionization scattering by creat- 
ing electron hole pairs in the first and the second bands 
respectively: eg = e2-\-e\+h\ and e\ = el+e2+h2, where 
the upper subscript Li is the angular momentum in units 
of 2/ 3d. In addition, angular momentum determines the 
band edge energy A^. In the large d limit, = hvpLi, 
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FIG. 1: (a) Impact ionization rates and (b) impact excitation 
rates for a (25,0) nanotube (d — 2.0 nm) as a function of en- 
ergy (measured from the bottom of the first conduction band 
Ai), for the first seven conduction bands shown in red, blue, 
green, cyan, magenta, light red, and black in increasing order 
of the band index. The vertical dashed lines correspond to 
the bottoms of the conduction bands A^ (i = 2, 3. ..7) relative 
to Ai. (c) The vertical lines are allowed k-points of zig-zag 
tube with mod{n, 3) = 1 for the first four doubly degenerate 
bands ef\ The angular momenta Li, in units of 2/3d, mea- 
sures a minimum distance from the K {K') point of graphene 
to the ID lines of allowed k-points of a nanotube. 



where vp is the Fermi velocity of graphene. Therefore, in 
this limit, the energy and momentum conservations are 
simultaneously satisfied for the electrons at the bottom 
of the third and the fourth bands. However, for finite 
diameter tubes: A3 < A2 4- 2Ai in mod{n — m, 3) = 1 
chirality and A4 < Ai -I- 2A2 in mod{n — m, 3) = —1 chi- 
rality. This pushes the lowest impact ionization threshold 
energy Eths of a (25, 0) CNT in the third band at energies 
higher than its bottom A3, i. e. Eth3 > A3, and a power 
law dependence is obtained Ta^ (x (E — Eths)"' with 
a « 2. Whereas at the bottom of the fourth band A4, 
energy and momentum conservation laws are simultane- 
ously satisfied such that the threshold energy EtM coin- 
cides with the bottom of the band A4, i. e. EthA — A4. 
The energy dependence of the scattering rate at the on- 
set in the fourth band is similar to that of a step function 
rii4 oc Q{E — Ethi)- In a (26, 0) tube the impact ion- 
ization threshold in the third band coincides with A3, i. 
e. Eth3 — A3, and it is higher than A4 in the fourth 
band, i. e. Ethi > A4. For impact ionization to take 
place in the first (second) band via the decay channel 
e\ = el+el + /ij"^ (e| = + e\ + h^'^), the threshold 
energy must be at least Ethi > 3Ai {Eth2 > A2 + 2Ai) 
and momentum conservation along the ID wavevector 
in non-parabolic bands brings the threshold for the first 
(second) bands to energies beyond the scale of Fig. ^ 

In impact excitation, the exciton state formed has an 
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FIG. 2: Electron distribution function in a (25,0) nanotube at 
applied electric field of (a) F = 0.1 MV/cm and (b) F = 0.05 
MV/cm as a function of energy in the first six conduction 
bands shown in red, blue, green, cyan, magenta, and light 
red, respectively. Black curves are the total distributions. 
In (a) solid curves use impact excitation and dotted curves 
impact ionization scattering rates with phonons at T = 300 
K. In (b) we use impact excitation scattering, solid curves 
use phonon scattering at T = 300 K and dotted curves use 
acoustic phonon scattering at Tac = 300 K and optical phonon 
scattering at Top — 1500 K. 

effective mass enhanced by 30% over the free particle 
value (me + ruh), due to the Coulomb interaction [isf. 
(We find that the magnitude of the mass enhancement 
depends on the strength of Couloumb interaction as 1/e, 
independent of tube chirality.) This mass enhancement 
lowers the threshold energy. Indeed, in the limit of the 
infinite effective mass, the exciton dispersion is flat, and 
therefore, energy and angular momentum conservations 
alone would determine the impact excitation threshold. 
Thus, the impact excitation thresholds of the first two 
conduction bands (Fig.Q^) are significantly reduced from 
the impact ionization values . At the onset of impact 
excitation, the third and the fourth band electrons can 
decay as e| = el + Efi or e| — e\ + Ef2 and el = e| + i?i2 
or 64 = e\+E22- The relative efRciencies of the two decay 
channels depend on tube chirality and electron energy. 

In electro-optical applications, we need to know the ef- 
fect of impact excitation scattering on the excited state 
production rate. This depends on the hot carrier dis- 
tribution function, which we calculate by solving the 
steady-state multi-band Boltzmann equation in the pres- 
ence of an electric field, electron-phonon scattering mod- 
eled as in [T^ , and either impact ionization Eq. or im- 
pact excitation scattering Eq. (0) . The results are shown 
in Fig. 12 At energies above the optical phonon energy of 
about 180 meV, the carrier distribution is determined by 
the field F and the optical phonon mean free path Aop. 
The probability that an electron is accelerated to an en- 
ergy E = eF\ over a length A without being scattered by 



FIG. 3: Exciton production rate P per unit carrier in a (25,0) 
tube as function of the electric field for acoustic phonon scat- 
tering at Tac = 300 K and different optical phonon tempera- 
tures (from bottom to top) Top = 300 K (red), 900 K (blue), 
1200 K (green), 1500 K (cyan), 1800 K (magenta), 2100 K 
(light red) along with the best fit to Eq. (|SJ| for Aop = 25 nm, 
Et = 0.57 eV and Po ~ 18 ps~^. The inset shows scaling of 
all the calculations with the effective temperature in accord 
with Eq. 0. 

an optical phonon is exp(— A/Aop). Therefore, the car- 
rier distribution g, for energies above the optical phonon 
and below the impact excitation (ionization) threshold, 
is expected to follow g{E) oc exp {—E/eFXop), as seen in 
Fig. 121 Above the impact excitation (ionization) thresh- 
old, the scattering is an order of magnitude larger and so 
is the slope of In {g{E)) versus E, as in Fig.|21 The higher 
threshold for impact ionization results in a "hotter" car- 
rier distribution than in the case of impact excitation 
(Fig. lit). 

The non-equilibrium phonon effect on the carrier dis- 
tribution is shown in Fig.[2l3- Although at optical phonon 
temperature Top = 1500 K the phonon scattering rate is 
a factor 1.7 larger than at room temperature, the slope 
of In {g{E)) versus E, for energies below the impact exci- 
tation threshold, is similar to the Top — 300 K case. On 
the other hand, the hot carrier tail above the impact ex- 
citation threshold significantly increases with Top, which 
in turn increases the exciton production rate. 

The exciton production rate itself is given by the prod- 
uct of the carrier distribution, as in Fig. [3 and the impact 
excitation scattering rate, as in Fig. and it is shown 
in Fig. U] We find that the rate P can be well fitted by 
the following equation: 

P = Poexp = I , (5) 

V ^ {ksTopf + {eFKpf ) 

where Et is the impact excitation threshold energy and 
Pq is a constant. The form of Eq. |SJ| resembles the Boltz- 
mann distribution with an effective temperature in the 
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FIG. 4: Energy distribution of excitons produced in the 
impact excitation process in (a) a (25,0) tube and (b) a (26,0) 
tube with phonons at T = 300 K and different applied electric 
fields F in MV/cm; 0.03 - red, 0.06 - blue, 0.09 - green, 0.12 - 
cyan, 0.15 - magenta, and normalized to the total production 
rate P ~ J luiduj. 



presence of the field. The inset of Fig. O shows an ex- 
cellent scaling of the exciton production rate with the 
effective temperature. The diameter dependence of the 
fit parameters in Eq. ((Sj for 1 nm < d < 2 nm is very sim- 
ple: Eth ~ 1.22 (eV nm)/d, Xop « 14 x d, while Pq « 18 
ps~^ is nearly diameter independent. 

Excitons in nanotubes can be created either by op- 
tical pumping as in photoluminescence experiments, or 
by electrical excitation as in electroluminescence. In the 
photoluminescence case, the energy of the created ex- 
citons follows the photon energy of the excitation light 
source, and the angular momentum distribution is de- 
termined by the dipole selection rules. Whereas in the 
impact excitation process, the energy distribution of the 
generated excitons depends on the applied bias and it is 
not subject to the same selection rules. In Fig. 01 we plot 
the distribution function of the initially created excitons 
at different fields. The distribution has several distinct 
peaks, which we identify by analyzing the contributions 
from excitons with different angular momenta. At the 
onset of impact excitation, the third and the fourth band 
electrons generate the low energy excitons with finite an- 
gular momenta: Ef^, E12, E22, whereas at higher bias 
electrons in the first band produce unbound excitons of 
zero angular momentum, which contribute significantly 
to a broad continuum at energy of about 0.9 eV in Fig.^ 

We conclude that CNTs and other ID systems are 



characterized by unusually high impact excitation rates 
and can be used as efficient exciton sources to be em- 
ployed in fundamental studies of interacting boson sys- 
tems at high densities and in opto-electronic device ap- 
plications. 
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